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Simulation of Wilson fermion at finite isospin density* 

A. Nakamura a and T. Takaishi b 

a RIISE, Hiroshima University, Higashi-Hiroshima 739-8521, Japan 
b Hiroshima University of Economics, Hiroshima 731-0192, Japan 

For QCD with Wilson fermions at isospin chemical potential we study the finite phase transition on an 8 3 x 4 
lattice at k = 0.15. We use two gauge actions: Wilson action and DBW2 action. Both actions give the same 
results. The phase diagram is qualitatively similar to the one obtained for QCD at small baryon chemical potential. 
We also calculate the number density for various isospin chemical potentials. 



1. INTRODUCTION 

Lattice simulations of finite density QCD are 
difficult due to the sign problem. Namely the 
action is complex and the complex phase fluctu- 
ates, which makes the simulations difficult. Re- 
cently considerable progress has been made for 
small baryon chemical potential /LtsPP- There ex- 
ist several approaches to small \±b '■ Taylor expan- 
sion method|2], reweighting method^ and imag- 
inary chemical potential method^]. 

Apart from the complex action, a model with 
a positive measure like isospin model can be used 
to obtain insights to QCD at finite hb- More- 
over it might be expected that the phase diagram 
of QCD at small /is is similar to that at small 
isospin chemical potential jii 5 . 

In this study we use Wilson fermions with fj,j 
and study the phase diagram and the number 
density. 

2. ISOSPIN DENSITY 

Lattice QCD partition function with Nf flavors 
is given by 

Z= [dU]{Y[detD(^))exp(S g [U}) (1) 

i 

where 

S g [U] = |((l-8ci) J2 ReTWIl+ci ReTrW12)(2) 



and D(fj,i) is the Wilson fermion matrix at /i^. 
Wll and TU12 stand for lxl and 1x2 loops 
respectively. For Nf = 2 with )i u = —fid = A*z we 
obtain 



det J D( y u u )detI)(^ d ) 



= det£>(/i/)detD(-/i/) 
= IdetDOij)! 2 (3) 



* Presented by T.Takaishi 



where a relation — 75-D ( — A*i)75 is used. 

We call fii isospin chemical potential. For this hi 
the measure is positive definite and the standard 
Monte Cairo technique can be applied. 

The parameter c\ in eq.@ is the one which 
specifies the gauge action. For example c\ = 
for the Wilson gauge action, and ci = —1.4089 
for the DBW2 action®. 



3. SIMULATIONS 

We generate configurations on an 8 3 x 4 lattice 
at k = 0.15 by the hybrid Monte Cairo (HMC) 
algorithm. We use the Wilson gauge action and 
the DBW2 action 6 . Simulation parameters are 
summarized in Tables 1 and 2. We choose a step 
size so that the acceptance of the HMC algorithm 
becomes 60 ~ 70% which gives the maximum per- 
formance of the HMC algorithm with the 2nd or- 
der leapfrog integrator^. This optimum accep- 
tance w 60 ~ 70% does not depend on the details 
of the action we take. It depends on the order of 
the leapfrog integrator. 
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Mi 


/? 


At 


Acc. 


Traj. 


0.0 


5.38 


1/14 


0.667(4) 


13000 


0.1 


5.37 


1/14 


0.667(5) 


12000 


0.2 


5.365 


1/14 


0.652(5) 


15500 


0.3 


5.34 


1/14 


0.647(5) 


12500 



Table 1 

Simulation parameters for the Wilson gauge ac- 
tion. At is the step size and the trajectory length 
is set to 1. "Acc." and "Traj." stand for Accep- 
tance and # of trajectories respectively. 



Mi 





At 


Acc. 


Traj. 


0.0 


0.68 


1/16 


0.644(3) 


28000 


0.1 


0.67 


1/16 


0.646(2) 


17800 


0.2 


0.66 


1/16 


0.653(3) 


30000 



Table 2 

Same as in Tabic 1 but for the DBW2 action. 



4. RESULTS 

In order to locate the phase transition point 
we measure susceptibility for various observables. 
Since for each \ii we performed a single simulation 
at one 0, the reweighting method|S] was used to 
investigate a region in the vicinity of the 0. The 
expectation value of an observable O at 0' can be 
obtained through a single simulation at by 



(0)p = 



(Oexp(A/?))g 



(4) 



where A0 = 0' - 0. 

Measurements are done for lxl, 1x2, 2x2 
Wilson and Polyakov loops. Figure ^ shows the 
lxl loop from the Wilson gauge action for various 
\ii . The critical coupling constant C decreases as 
\±i increases. The precise value of C is estimated 
by measuring susceptibilities of the observables. 

Typical examples of susceptibilities from vari- 
ous observables are shown in Figures EE All sus- 
ceptibilities from different observables give simi- 
lar C . From C (Wilson gauge action) evaluated 
from the lxl loop susceptibility we obtain 



c (hj) = 5.38 - 0.54^f . 



(5) 



The coefficient of O(fij) term is similar to those 
of 0|. Figure 01 shows the phase diagram in the 
T — fii plane. is converted to T by the 2 loop 
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Figure 1. lxl Wilson loop from the Wilson 
gauge action for various fii as a function of 0. 
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Figure 2. Susceptibility from various observables 
at fii = 0.3. The Wilson gauge action is used. 



function. Two phase diagrams from the Wilson 
and DBW2 gauge actions are in good agreement. 
The phase diagrams are qualitatively similar to 
the one obtained with KS fermions at small fis 3 

m. 

We also calculate the number density rid de- 
fined by — Tr — — — . Figure [5] shows rid, f° r dif- 
V D dfxi 

ferent fij. rid seems to increases with 0. 



5. DISCUSSION 

We have studied QCD at finite isospin density 
with Wilson fermions. The results show that the 
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Figure 3. Same as in Figure El but from the 
DBW2 action at ^ = 0.0. 
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Figure 5. Number density for various /ij as a 
function of (3. We generated 25—50 configurations 
at each (3. 
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Figure 4. Phase diagram in the T — \ij plane. T 
is normalized by Tc at /i/ = 0.0. 



phase diagram is qualitatively similar to the one 
obtained at small [ib- The Wilson and DBW2 
gauge actions give the same phase diagram. More 
quantitative analysis is needed to confirm the 
present results. 
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